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Abstract 

We analyze high-temperature series expansions of the two-point and four-point correlation- 
functions in the three-dimensional euclidean lattice scalar field theory with quartic self-coupling, 
which have been recently extended through twenty-fifth order for the simple-cubic and body- 
centered-cubic lattices. We conclude that the length of the present series is sufficient for a fairly 
accurate description of the critical behavior of the model and confirm the validity of universal- 
ity, scaling and hyperscaling. In the case of the body-centered-cubic lattice, we determine the 
value of the quartic self-coupling for which the leading corrections to scaling approximately vanish 
and correspondingly the universal critical parameters can be determined with high accuracy. In 
particular, for the susceptibility and the correlation-length exponents we find 7 = 1.2373(2) and 
v = 0.6301(2). For the four-point renormalized coupling we find g = 23.56(3). In the case of the 
simple-cubic lattice our results are consistent with earlier estimates. 
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I. INTRODUCTION 



High-temperature (HT) expansions for scalar spin models with bilinear nearest-neighbor 
interaction and general single-spin measure, on 2- and 3-dimensional bipartite lattices, have 
been extended ^ 2 i 3 i 4 i 5 i 6 in the last few decades from order K 10 to K 25 . Here K = denotes 
the usual HT expansion variable, with k the Boltzmann constant, T the temperature and 
J some energy scale characterizing the next-neighbor interaction. The main observables 
expanded, include the susceptibility X2, the second moment of the correlation-function /i 2 , 
the zero- momentum four-point function xa- The linked-cluster method used in these calcu- 
lations expresses the series coefficients in a closed form as polynomials in the (normalized) 
moments of the single-spin measure thus making the analysis of a broad class of models 
possible. This fact entails two related benefits. First, it enables one to explore the extent 
of the Ising universality class to which these models should generally belong. Second, one 
can take advantage of the non-universality of the amplitudes of the corrections to scaling in 
the critical behavior to obtain very accurate estimates of universal critical parameters, such 
as exponents and scaling functions, by focusing the numerical analysis on those particular 
spin models, within the given universality class, which show vanishing (or quite small) lead- 
ing corrections to scaling(LCS). This prescription, which can be effective only with rather 
long series, was suggested in Refs.-^ and extensively tested by 21-term HT series on the 
body-centered-cubic lattice, for various single-spin measures. 

In the initial studies, however, the available series were simply too short. In Ref.— an 
analysis of the critical properties of the lattice 4 euclidean field theory in two-, three- 
and four- dimensional space, using tenth-order HT series, suggested a complicated and puz- 
zling critical behavior in the 3d case, which showed failures of universality and hyperscaling 
and cast doubt on the validity of key mathematical assumptions^ ( see however Ref.—) in 
the application of the renormalization group to the critical phenomena. This analysis was 
historically important in raising questions which could be answered only by numerical inves- 
tigation of specific models and was a strong incentive to further extend the HT series and 
to repeat and improve simulation studies. The anomalous features observed in Ref.— were 
not confirmed by the successive studies^&i 1 . and were finally ascribed to the shortness of the 
series analyzed. Therefore over a decade ago the problem was considered as settled. 

We have recently extended, through order K 25 , the HT expansions of the moments of 
the correlation-function and, through order K 23 , the expansion of the zero- momentum four- 
point correlation-function for a general scalar model defined on three-dimensional bipartite 
lattices (such as the simple-cubic (sc) and the body-centered-cubic (bcc) lattices). These 
expansions have already proved useful for accurately displaying the properties of critical 
scaling and universality of the spin-S* Ising models**^ and have also produced high-precision 
estimates of critical exponents and universal combinations of critical amplitudes, based on 
the prescriptions of Refs.i^. Also following the lead of the same Authors^, analogous 
expansions, independently obtained^ for the sc lattice case, through order 25 for \2 and fi 2 
and through order 21 for x 4 , were used to study the lattice <p\ model in a neighborhood of 
a particular value of the self-coupling for which the LCS are negligible. It remains to show 
that the present HT expansions are also adequate for an extensive and accurate study of 
the scalar field model for two different lattices and on the whole range of the self-coupling. 
From our study we conclude that the present HT series are sufficiently long that, by properly 
resumming them, our analysis can reach an accuracy comparable or higher than the most 
extensive MonteCarlo simulations. Moreover, we show that also in the case of the bcc 
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lattice, we can approximately determine a value of the quartic self-coupling such that the 
LCS vanish, and thereby we can obtain^ very accurate estimates of the universal critical 
parameters in complete agreement with the other cited recent computations^. 

In the following paragraphs we shall first recall the definitions of the quantities whose 
HT expansions are studied. We shall then very briefly comment on the results of the series 
analysis. We shall map out the phase diagram of the model and exhibit its universality 
properties along the critical line, such as the independence of the critical exponents 7 and 
v on the self-coupling of the field and on the lattice structure. Then, we shall give evidence 
that the critical renormalized coupling constant is universal and non-zero, thus verifying the 
validity of hyperscaling. Finally we shall report our best estimates of the critical exponents 
7 and v and of the renormalized four-point coupling from the analysis of the model with 
minimal LCS on the bcc lattice. The figures summarizing our results are among the main 
motivations of this study, since some of them qualitatively differ from the analogous ones 
appearing in Ref.— , which, to the best of our knowledge, have never been revised in the 
literature. 



II. THE MODEL 

The lattice $\ model is defined by the Hamiltonian 

KH{<p} = -K £ Mj + + goiti - i) 2 ) (i) 

<i,j> i 

Here i and j are integer-component (multi) indices denoting the sites of a three-dimensional 
lattice and the first sum extends to all nearest-neighbor sites. 

This Hamiltonian is obtained from the lattice discretization of the euclidean action of the 
continuum real field ip(x) in 2 + 1 dimensions. 

S = J[\{d^f + \m^ + ^}dx (2) 

after setting ijj = m 2 = _ 6 and \ = <|a 

We will analyze the HT expansion of the moments of the two-point connected correlation- 
function < > c and of the zero-momentum connected four-point correlation-function 

X4,(K>9o) = T,i,j,k < 4>o<Pi4>j(pk >c- 

The susceptibility is defined by the zero-order moment of < 4>o4>i >c 

X2(K,g ) = £ < o 0i >c ■ (3) 

i 

Its behavior, as K tends from below to the critical value K c (go), is expected to be 

X 2(K,g ) = A(g )r(g )-^[1 + a(g )r e ^ + ...] (4) 

where A(g ) is the critical amplitude of the susceptibility, r(g Q ) = 1 — K/K c (g ) is called 
the reduced (inverse) temperature, 7(^0) is the critical exponent of the susceptibility, a(g ) 
is the amplitude of the LCS and 9 (go) is the exponent of the LCS. 

The square of the correlation-length is expressed in terms of the ratio of the second 
moment of the correlation-function ^(K, go) =J2i^ 2 < 0o0« >c and the susceptibility as 

? \K, 90) = 7 777 7 ( 5 ) 

6x21^, go) 
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and is expected to show the critical behavior 



e(K,g ) = B{g G )r{g,)- 2 ^[l + b(g o )r ^ + ...] (6) 

where v(go) is the critical exponent of the correlation-length. 

The critical renormalized coupling constant g(go) is expressed in terms of Xi(K,go), 
X2(K, g ) and ^ 2 (K, g ) as the value of the quantity 

(tr \ ~ v Xa(K, g ) 

when K tends from below to K c (go). Here v denotes the volume per lattice site. We have 
v — 1 for the sc lattice and v = 4/3 v^3 for the bcc lattice. The expected critical behavior of 
g(K,g ) is 

g(K, g ) = ^(^o)r(^o) 7(9o)+3l/(9o) - 2A4(9o) [l + c(g )r e ^ + ...] (8) 

where A 4 (g ) is the gap exponent. If j(go) + 3v(g ) — 2A 4 (g ) = 0, we say that hyperscal- 
ing holds and the critical renormalized coupling g(go) is the finite positive 13 critical limit 
of g(K,g ). Of course, the exponents 7(^0), v(g ), A 4 ((jo), 9 (go) an d the critical renormal- 
ized coupling g(go), as determined from the analysis of the HT series, will appear to be 
independent of g , if universality is valid along the critical line. 



III. NUMERICAL RESULTS 
A. The phase diagrams 

In Figs. [U and El we have mapped out the phase diagrams of the 0| model on the sc 
and the bcc lattices, respectively. There are two phases separated by a line K = K c (g ) of 
second-order critical points. In the disordered phase, below the critical line, the reflection 
symmetry <p — ► —<p is unbroken. For g = the model reduces to the Gaussian model and 
we have K c (0) = 2/q = 1/3 on the sc lattice (K c (0) = 2/q = 1/4 on the bcc lattice), where 
q is the lattice coordination number. For g — > 00 the model reduces to the usual spin-1/2 
Ising model and we have K c (oo) = 0.221655(2) on the sc lattice (K c (oo) = 0.1573725(10) 
on the bcc lattice) as indicated in our general study 5 of the spin-S 1 Ising model. In Table 
1 we have reported numerical estimates of K c (g ) for a few values of < g < 00 in the sc 
and the bcc lattice cases. These estimates are obtained by the following simple procedure. 
Using the susceptibility expansion, for each value of go, we form the sequence {K^ n \g )} of 
Zinn- Justin modified-ratio approximants (MRA)2*i^ of the critical point. We observe that 
these sequences approach smoothly their expected asymptotic behavior— 

Kjr\g ) = K c (g )[l - C( 7 , 9)a(g )/n 1+e ^ + o(l/n 1+e ^)] (9) 

where a(g ) and 9 (go) are, respectively, the amplitude and the exponent of the LCS appearing 
in eq.(JH) and C(j, 9) is a positive constant 5 - depending on the values of 7(^0) an d 9(g ). We 
expect that the exponent 9(g ) be universal, namely independent of g and indeed, to a 
good approximation, the asymptotic behavior of the MRA sequences is consistent with the 
value 9 = 0.517(4) suggested by a recent simultaneous studyi 5 - of a set of models in the 
Ising universality class. Since, for all values of go, the contributions of the higher-order 
terms in eq.© appear to be small, particularly so in the bcc lattice case, we shall assume 
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that the MRA sequences {K^(g )} can be extrapolated to infinite length n of the series 
simply by fitting the highest-order (alternate) approximants to the first two terms of the 
asymptotic expansion (JOJ). The values of K c (g ) determined by this prescription are quite 
stable, both under small variations of 9 and of the fitting procedure, such as including in 
the fit higher-order terms of the expansion Q. Our conclusion is that at least five decimal 
figures of the result are reliable. These estimates are then refined by comparing the results 
of the extrapolations with the determination of K c (g ) from first-, second- and third-order 
differential approximants(DA)^*^ and finally error bars are inferred which are likely to be 
only upper bounds. As expected, these uncertainties, which reflect the differences between 
the estimates by extrapolated MRA sequences and by DAs, depend strongly on the size of 
the corrections to scaling and therefore on the value of go, but of course they are generally 
invisible on the scale of figures and EJ 

B. The critical exponents 

In Figs. El and El referring to the sc and the bcc lattices respectively, we have plotted our 
estimates of the exponent of the susceptibility 7(^0) vs - Gl = Also in this case, using 
the susceptibility HT expansion, for each value of g , we can form the sequences {7^ (#0)} 
of MRAs of the exponent 7 (go) an d compare them to their expected asymptotic expansion 5 

7 (n) G?o) = l(9o) - D( % 9)a(g )/n e ^) + (l/ n «C»)). (10) 

Here D(j, 6) is a positive constant^ depending on 7(^0) an d 9 (go) and a(g ) is the amplitude 
of the LCS in eq.(jU). We can observe that also these MRA sequences are smooth, but a simple 
minded extrapolation based on a fit of the sequences {7^ (go)} to the first two terms of the 
asymptotic expansion (|TUJ) cannot lead to very accurate estimates of j(go), because, even at 
the present order of expansion, the contribution of higher-order corrections is not sufficiently 
small. In the figures El and EJ we have represented by solid lines the estimates 7 (go) obtained 
from the MRAs 7 (l7 ' ) (go), I^Hgo) and 7 < - 25 - ) (go), which use the susceptibility series only up 
to the orders 17, 21 and 25 respectively. These estimates show a rapid crossover from the 
Gaussian value 7(0) = 1 to a behavior which, for g > 0, tends to become independent of go 
and to approach the dashed line in the figure indicating the central value 7 = 1.2371 of our 
estimate in Ref.— for the Ising universality class. If we include higher correction terms in 
the asymptotic expansion (fTU|) somewhat more accurate estimates of 7 (go) can be obtained. 
The simplest possibility is to introduce a single effective higher correction 0(1/ n s ) with 
go-independent exponent s. Choosing the value s = 4, we obtain the estimates represented 
by the dotted curve. The accuracy of the approximation can be further improved at the 
expense of the simplicity of the fitting procedure, but here it is sufficient to indicate only 
the qualitative trend. In the same figures, we have also reported estimates of 7(^0) obtained 
by simple first-order unbiased DAs using all available series coefficients. In order to keep 
the figures readable, these data are restricted to a smaller range of values of g of particular 
interest. Completely analogous results are obtained starting with the HT series for the 
correlation- length squared £ 2 to determine the critical exponent v but, for brevity, they will 
not be reported here. 
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C. Leading corrections to scaling and accurate estimates 

In Fig. El we have plotted vs. go the amplitude a(go) of the LCS in as obtained by 
fitting the MRA sequence {K^(g )} to the asymptotic expansion Q in the bcc lattice case. 
Of course, the determination of a(go) is much more sensitive than K c (g ) to the inclusion 
of higher-order corrections in the fit to the expansion (j5J). The most convincing results are 
obtained when a single additional correction 0(1/ n e+0 ) is included. This three-parameter 
fit changes the results of the two-parameter fit for K c (go) only within the estimated uncer- 
tainties and does not alter the qualitative structure of a(go), but only brings its evaluation 
into closer agreement with the corresponding estimate of a(go) from a similarly improved 
fit to eq. ()10p and with the values of a(oo) obtained in earlier works^. The uncertainty of 
the final results should not exceed 5 — 10%. In the case of the sc lattice, the behavior of 
a(go) is similar, but a(g ) cannot be determined with comparable accuracy by this straight- 
forward method because it is even more sensitive to the presence of higher-order corrections 
in 0. The important fact is however that, both in the sc and the bcc lattice cases, when 
g increases from to oo, the amplitude a(g ) of the LCS in eq.(JU) varies from positive to 
negative values. In the sc lattice case a(go) vanishes at §q c ~ 1.10(2). We recall that the 
above cited numerical analysis of the 0| HT series on the sc lattice^ was indeed performed 
precisely at go = 1.10. In the bcc lattice case we find that the zero of a(go) occurs at 
g^ cc pd 1.85(5). In Ref.— cjq c had been determined by an appropriate MonteCarlo simulation. 
In our approach no preliminary MonteCarlo simulation is needed to determine g^ 00 , since, 
knowing the structure of the function a(g ), our study of the model on the whole range of 
values of g directly gives a sufficiently accurate indication. For g in a neighborhood of the 
zero of the LCS amplitude, we can observe that the various approximations of 7(^0)? shown 
in Figs. E3 and EJ tend to coincide and come nearest to our earlier estimate^ for the spin S 
Ising model: 7 = 1.2371(1) shown in the figure by a dashed line. Eq. fjlOj) and the measured 
behavior of a(go) simply explain why, for g < c/q c , (or g < g h Q cc in the bcc lattice case) the 
MRAs of Figs. H3 and 0] approach their constant limiting value from below and otherwise 
from above. Of course, the same mode of approach to the limit can be observed for the 
estimates obtained from DAs, when varying the number of HT series coefficients used in the 
approximants, however, for clarity we have reported in the figures only the highest order 
results. 

We can now take advantage of our estimate of g h Q cc to compute some universal critical 
parameters of the (f)\ model on the bcc lattice with approximately vanishing LCS. The HT 
series coefficients of X2, ^2 and %4 011 the bcc lattice, for g = 9o cc , are reported in Table 
II. (The corresponding expansions in the sc lattice case, for go = g'o c , can be found in 
Ref.—.) Very accurate estimates are thus obtained for the critical exponents and the critical 
renormalized coupling. In particular, using first- and second-order DAs (either unbiased or 
biased with the critical value K c (g b cc ) = 0.2441357(5)), we find 7 = 1.2373(2), v = 0.6301(2), 
g = 23.56(3) in good agreement with our previous estimates^ and with the results of Ref.—. 
Our estimated errors account also for the ~ 3% uncertainty in the estimate of g^ 00 . 

Since our HT series for g(K, go) is two terms longer, we have also repeated the analysis of 
Ref.— with our sc lattice series computed at ^ c = 110. The final estimates: 7 = 1.2372(2), 
v = 0.6301(2) and g = 23.56(4), evaluated at K c (g s c ) = 0.375097(1), are consistent both 
with our results for the bcc lattice and with those of Ref.—. 
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D. The renormalized coupling 



In Figs. Eland[7[ which refer to the sc and the bcc lattices respectively, we have plotted 
the renormalized four-point coupling constant <7(£ 2 ,<7o) vs - A 4 = -gt+i for several fixed 
values of g . For fixed g and £ 2 — > oo, all curves appear to tend to a g - and lattice- 
independent limiting value g(oo, go) which is consistent with our estimate in Ref.— of the 
critical renormalized coupling in the Ising universality class g = 23.52(5) indicated in the 
figure by the dashed line. This fact suggests that the critical renormalized coupling is 
universal with respect to the self-coupling go and to the structure of the lattice and moreover 
that it is non-zero (so that hyperscaling is valid). The curves in the figures 6 and 7 are 
obtained in a most straightforward way from the highest-order non-defective "simplified 
differential approximants" (SDA)^ biased with the values of 9 and K c (g ). Completely 
equivalent results are obtained also, slightly more laboriously, using ordinary biased DAs 
instead of SDAs. The estimated width of the error bars is comparable to the thickness of the 
lines. Also for g(£, 2 ,go), we can observe that, as g varies from to oo, the amplitude c(g ) 
of the LCS in (jSJ), which is negative for small g , changes sign at g , as it should, because 
the ratio c(go)/a(go) is expected to be universal. Correspondingly, as shown clearly by these 
figures, the curves tend to their common limiting value from below if go < go and otherwise 
from above. 

Figs. |H1 and 03 for the sc and the bcc lattices respectively, offer a different view of the 
same data. Here the renormalized coupling g(£ 2 ,go) is plotted vs. Gl = ^-j- for several 
fixed values of £ 2 . All curves are monotonically increasing in go and show a rapid crossover 
from the Gaussian value g(£ 2 , 0) = to a shape increasingly flatter as £ 2 becomes large. 
For g < %, as £ 2 — > oo, the sequence of curves tends from below to the constant g(oo, g ), 
for g > go the same value is approached from above. We have already stressed, that this 
pattern of behavior simply reflects the fact that the amplitude c(go) of the LCS in g(K,go) 
is negative for g Q < g and positive for go > go- Here we should also recall that in a similar 
plot reported in Ref.— (only the case of the bcc lattice case is discussed in that study) the 
curves show quite a different structure for £ 2 > 64. In particular, the curves g(£, 2 ,go) in 
Ref.— are not monotonic in g and their limiting behavior as £ 2 — > oo does not appear to be 
universal and non-zero. We can see no reason for these anomalous features other than the 
insufficient length of the ten-term HT series used in the analysis of Ref.—. 

One more view of the same data is presented in Figs. El and ^2 f° r the sc and the bcc 
lattices respectively, showing contour plots of g(£ 2 ,<?o) in the £ 2 = go = plane. 

Also these figures are qualitatively different from the corresponding ones of Ref.— which show 
a spurious saddle-point structure and suggest a failure of universality. Nothing like that can 
be inferred from our updated figures. 

E. Conclusions 

In conclusion, we have analyzed HT expansions for a one-parameter family of continuous 
spin models which interpolate between the Gaussian and the spin-1/2 Ising models. In 
the case of the bcc lattice, we have taken advantage of the parameter dependence of the 
amplitudes of the LCS to improve the accuracy in the determination of some universal 
critical quantities. Moreover, we have shown that, in the light of our extended series, all 
puzzling and unexpected features which emerged from the old HT analysis of Ref.— can only 
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be ascribed to numerical inaccuracies deriving from the use of too short series. 
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FIG. 3: The susceptibility exponent 7(50) as a function of Gl = for the sc lattice. The 
continuous curves represent the estimates obtained by the MRAs of orders 17, 21 and 25 respec- 
tively. The dashed curve indicates the central value of our estimate 5 of 7 for the Ising universality 
class. The dotted curve shows the results of an extrapolation of the MRA sequences by a fit to 
the asymptotic expansion H10|) supplemented by a single higher-order correction 0(l/n 4 ). The 
triangles show estimates of 7(170) by simple unbiased first-order DAs. 





FIG. 5: The amplitude a (go) of the leading correction to scaling in (J2J) vs. Gl = , as 
obtained from a fit to the asymptotic expansion @ supplemented by a single higher-order correction 
0(l/n 6+e ) in the bcc lattice case. 




[ i i i i i i i i i i i i i i i i i i i ] 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 6: The renormalized coupling g(£ 2 ,go) as function of XA = for several fixed values of 
the (j) 4 self-coupling go on the sc lattice. The values of go increase from the lowest curve up. The 
dashed line represents the central value of our previous estimate g = 23.52(5), obtained in Ref£, 
of the renormalized coupling in the Ising universality class. 
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I i □ r~ i i i i i i i i i i i i i i i i 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 7: The same as in Fig|f2 but for the bcc lattice. 




FIG. 8: The renormalized coupling g(£, 2 ,go) as function of Gl = for several values of the 
correlation-length square £ 2 on the sc lattice. The values of £ 2 are indicated by the symbols 
attached to the corresponding curves. The asterisk labels the curve for £ 2 = 1, the empty rhomb 
for £ 2 = 2, the full rhomb for £ 2 = 4, the empty circle for £ 2 = 9, the full circle for £ 2 = 16, the 
empty square for £ 2 = 36, the full square for £ 2 = 100, the empty triangle for £ 2 = 256, the full 
triangle for £ 2 = 900. The dashed line represents the central value of our estimate g = 23.52(5), 
obtained in Ref. 5 of the renormalized coupling in the Ising universality class. 
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0.0 0.2 0.4 0.6 0.8 1.0 

G1 

FIG. 9: The same as in FigEl but for the bcc lattice. 




FIG. 10: Contour plot of the renormalized coupling g = g((, 2 ,go) in the G4, X4 plane for the sc 
lattice. Here G4 = anc j X4 = |iqqt • The values of g{£, 2 ,go) are shown beside the corresponding 
level curves. 
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FIG. 11: The same as in Fig llUl but for the bcc lattice. 
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TABLE I: K c (go) for a few values of of go in the case of the sc and bcc lattices. 



5o 



K s c c (g ) 



K 



bcc I 



0.000 0.3333333 

0.010 0.34055(1) 

0.050 0.35940(1) 

0.100 0.37340(1) 

0.150 0.38238(1) 

0.200 0.38842(1) 

0.250 0.39250(1) 

0.300 0.39517(1) 

0.350 0.396782(4 

0.400 0.397585(3 

0.500 0.397397(3 

0.650 0.394135(3 

0.760 0.390343(3 

0.900 0.384503(3 

1.060 0.377035(2 

1.100 0.375096(2 

1.140 0.373137(3 

1.300 0.365220(3 

1.600 0.350584(3 

1.800 0.341319(3 

1.850 0.339085(3 

2.000 0.332596(3 

3.000 0.297990(3 

4.100 0.274459(3 

5.000 0.262538(3 

6.400 0.251303(3 
10.000 0.238853(4 
15.000 0.232587(4 
45.000 0.225112(4 



0.250000 
0.25471(1) 
0.26724(1) 
0.27655(1) 
0.28242(1) 
0.28627(1) 
0.28876(1) 
0.29029(1) 
0.29108(1) 
0.29131(1) 
0.29055(1) 
0.28739(1) 
0.284131(4) 
0.279335(3) 
0.273375(3) 
0.271844(2) 
0.270306(2) 
0.264137(2) 
0.252883(1) 
0.245830(1) 
0.2441357(5) 
0.239230(1) 
0.213328(1) 
0.195919(1) 
0.187155(2) 
0.178926(2) 
0.169848(2) 
0.165296(2) 
0.159875(2) 
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TABLE II: The HT expansions of the susceptibility X2> the second moment of the correlation- 
function [j,2 and the zero-momentum four-point function X4 evaluated at go = 1.85 on the bcc 

lattice. 

order %2 ^2 Xa 

0.607989217187788034965 0.000000000000000000000 -0.530149722377688222012 

1 2.957207105732954317729 2.957207105732954317729 -10.31441007042348207915 

2 13.09429900741101048493 28.76720053242789148965 -116.8619663540111359343 

3 57.79323569148811046867 197.7144175506477337902 -1037.709279795438675000 

4 248.2179229506115615651 1185.408945561549742349 -7930.113134914392838073 

5 1064.858590911481926754 6551.599918767001108148 -54882.22830020753646424 

6 4512.534104619384245591 34421.02131513359339522 -353349.3513943436529616 

7 19111.25951155995320831 174321.7172035645509269 -2155038.392864457644076 

8 80383.21948173478800467 859681.0962551252453637 -12593879.17023276266910 

9 337958.2776471983168127 4151521.352493417270411 -71122511.16749072143498 

10 1414677.216907545784892 19721977.30824040566487 -390463929.5724427479526 

11 5919934.511772996630729 92423661.55059636566134 -2093673545.070737911244 

12 24698142.58226079299701 428322758.4354179679972 -11002848880.38386524575 

13 103016201.6436073223436 1966128942.161860887810 -56833474652.32522843381 

14 428735652.9640138248084 8952466557.867528601072 -289181057575.1432464575 

15 1783970290.009394230477 40476319807.17866188775 -1452144744869.861166239 

16 7410672453.259276817146 181885433500.9131802296 -7207354502435.721114810 

17 30779127150.29510373162 812879825813.3598744239 -35401801340419.05097854 

18 127668343287.7040888522 3615487616835.276408009 -172273869711328.5865152 

19 529480497190.9511942304 16011252152185.13150869 -831302587175201.5552329 

20 2193593476094.224975766 70632026759664.65740219 -3980914502340770.007282 

21 9086796325630.050694943 310488939709344.3260940 -18931663273718966.90898 

22 37608562966496.06439966 1360524482339925.373685 -89460744874166538.46043 

23 155638849919474.9034983 5944225695917982.872040 -420283001120399961.6139 

24 643624490532272.2532159 25901506033628787.27895 

25 2661386657660504.869714 112585550276346019.8717 
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